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The muon-catalyzed-fusion processes in D2 gas for various temperatures and den- 

r^ [ sities have been studied by means of Monte Carlo simulations. In particular, the role 

of the resonant and nonresonant ddfi formation and differences between the neutron 

time spectra from the ortho-D2 and para-D2 targets at low temperatures have been 

investigated. 
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Abstract 



• ,-{ • The negative muons which slow down in a D2 target form the dfx muonic hydrogen 

j>.^. atoms in the highly excited Rydberg states. Then, such atoms undergoing the various 

^ ! cascade processes [I] deexcite to the ground IS state and participate in numerous processes 

■^' accompanying fiCF (muon-catalyzed-fusion) phenomenon [2]. The ground-state IS atoms 

can gain high kinetic energies up to ~ 1000 eV, as a result of several accelerating collision 

>■ ■ processes during the atomic cascade (see, e.g., Refs. [3llll[5l[6]). However, a fraction of the 

^ ■ IS muonic atoms is thermalized after the cascade. 

The aim of this work is to study the role of "hot" dfi atoms, supplied by the cascade, 
and the competition between the energy-dependent processes (muonic atom scattering and 
muonic molecule formation) of /iCF cycle in a pure deuterium D2 gas. 

The D2 targets simulated in our Monte Carlo calculations correspond to the D2 targets 
O ■ which are to be used at JINR in the upcoming experiment [7j on the radiative capture 
reaction d + d — »■ ^He + 7 in the ddfi muonic molecule. The deceleration, thermalization and 
►v> I diffusion of d/i's in D2 perfect gas is described using the differential scattering cross sections. 
^ I At collision energies e greater than about 1 eV, the doubled cross section for dfi scattering 
from the bare deuterium nuclei are employed [H [9]. At lower energies, where molecular 
binding and electron-screening effects are significant, the differential cross sections for dfi 
scattering from isolated D2 molecules are used [10]. The total "molecular" and "nuclear" 
cross sections are smoothly sewed together, which was presented in Ref. jTT]. The cross 
sections used take into account the spin-flip of dfi atom and the rotational- vibrational 
excitations of the target molecules. For the considered temperatures T = 40-300 K, the 
D2 molecules are in the ground vibrational state z/ = 1 . The population of the total spin 
states F = 1/2 and 3/2 after the cascade is assumed to be statistical. 

The dfi atoms can form the muonic molecules ddfi in collision with the D2 molecules. 
Such a process can have both the nonresonant or resonant character [2]. The latter is 
possible due to the presence of the loosely-bound state rotational-vibrational state J = 



f = 1 of ddfi Ll2j. The energy excess in nonresonant dd/i formation is taken away by the 
electron [13] 

(i/i + D2 — ^ [{dd^)de]'^ + e . 

The rate of this reaction is practically independent of the dfi atom spin. At e < 100 eV, 
only the s and p waves give significant contributions to the formation rate. In Fig. [H the 
calculated rate of nonresonant dd^ formation in dfi collision with a free D2 molecule is 
plotted versus d/i energy [13]. The formation rates in this paper are given for the target 
density = 1 in the liquid hydrogen density units (LHD). The nonresonant ddfi formation 
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Figure 1: The rate of nonresonant ddfi formation in dfi scattering from D2 gas versus dfi 
energy. 

rate is significant (~ l/is~^) even at a few tens eV, has maximum at £ ~ 7 eV, and falls 
rapidly when e ^ 0. 
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Figure 2: The effective rate of resonant dd/j, formation in dfi scattering from 40-K ortho-D2 
{K = 0) and para-D2 {K = 1) gas versus dfi energy. The solid lines correspond to the 
upper spin state F = 3/2 of dfi, the dashed lines describe the F = 1/2 state. 

In the resonant ddfi formation, the energy excess is transferred to rotational-vibrational 
excitations of the muonic-molecular complex 



dfi + D2 ^ [{ddfi)dee]* . 



The resonant rate strongly depends on the d^ energy and its total spin state F. In Fig. [21 
the effective resonant ddfj, rate (with the back decay of dd/j, included) for the perfect D2-gas 
target at 40 K, calculated according to the method presented in Ref. [Hj, is shown for the 
states F = 3/2 and 1/2. Comparing Figs. [T] and [21 one can see that the resonant formation 
dominates at lowest energies (e < 0.1 eV). The resonant rates calculated for the rotational 
states K = and 1 of the target D2 molecule display significant differences. As a result, 
the time spectra of d-d fusion products should expect appreciable ortho-para effect at low 
temperatures. 

The Monte Carlo simulations for the 40-K D2 target have been performed using the 
nonresonant ddfi formation rates shown in Fig. [1] The resonant dd/i formation has been 
described using the absolute resonant formation rates and the back-decay rates. Since the 
mean free path of (i/i's is much smaller than the dimensions of the JINR target (the volume 
of 180 cm^), the simulations have been carried out for the infinite D2 target. The nitrogen 
contamination of 10~^ have been taken into account. The reliable theoretical distribution 
of dfi energy after the cascade process is not calculated yet. Therefore, following Ref. [3], 
the initial d/i energy is assumed to be a sum of the two Maxwell distributions. The first 
Maxwell component corresponds to the nonthermalized rf/x's and the second one describes 
the thermalized atoms. A relative population of these two groups of rfyu's is used as a free 
parameter in our simulations. 

In Figs. [3H6] are shown some results of our calculations for the 40-K D2 gas. The mean 
energy £avg of dfi's in the both spin states is shown in Fig. [3] for short times. Only the 
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Figure 3: Mean kinetic energy of d^{F = 1/2) and dfi{F = 3/2) atoms in D2 gas (T = 40 K, 
density = 0.05) versus time. The initial dfi energy distribution has the Maxwellian form 
with the mean energy of 10 eV. 

nonthermalized Maxwell component with e^vg = 10 eV has been considered, in order to 
study the dfi deceleration from high energies. One sees that slowing down to 0.05 eV is 
very fast — it takes about 25 ns. Then, this process is much slower since dfi's can gain 
appreciable energy due to the thermal motion of the target D2 molecules. At this stage 
of thermalization a difference between the mean energies in the states F = 1/2 and 3/2 is 
apparent. This is caused by the downwards spin-flip reaction 



dfi{F = 3/2) + D2 ^ rf/i(F = 1/2) + D2 



in which the energy Ai^hfs = 0.0495 eV (in the dfi + d center of mass) is released. The 
dash-dotted hne shows the thermal energy of 0.005 eV corresponding to T = 40 K. The 
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Figure 4: Spectrum of neutrons from d-d fusion in D2 gas (T = 40 K, density (p = 0.05) for 
short times. The contributions from all ddjj, formation processes and from resonant process 
are shown separately. 



time spectra of neutrons from muon-catalyzed d-d fusion for the same initial conditions 
are shown in Fig. |H The solid line represents the time spectrum with both the resonant 
and nonresonant ddfi formation processes taken into account. The dashed line has been 
obtained assuming only the presence of the resonant formation. At short times (< 20 ns) 
the fusion is mainly due to the nonresonant ddfi formation since most c?/i's are not slowed 
down to the energies corresponding to the resonance peaks. This phenomenon is more 
pronounced at lower target densities 
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Figure 5: Mean kinetic energy of dn{F = 1/2) and dn{F = 3/2) atoms in ortho-D2 and 
para-D2 gas (T = 40 K, density (j) = 0.2) versus time. The initial dfi energy distribution 
has the Maxwellian form with the mean energy of 10 eV. 

When the temperature of a 300-K D2 gas is lowered, the even rotational states deexcite 
to i^ = and the odd rotational states fall to K = 1, because of the symmetry of the 



D2 wave function. As a result, the cooled deuterium target forms the statistical mixture 
(2:1) of the ortho-D2 and para-D2 molecules [TB]. The deceleration of dfi{F = 3/2) and 
diJ,{F = 1/2) atoms in ortho-D2 and para-D2 targets at T = 40 K and = 0.2 is illustrated 
in Fig. O The ground-state dfi{F = 1/2) atoms are fully thermalized after about 1000 ns. 
On the other hand, the mean energy of the dfi{F = 3/2) atoms never reaches the thermal 
value. The initial population of these excited atoms disappears after a few hundreds ns. 
For larger times, one sees only a very small amount of the dfi{F = 3/2) atoms that are 
created via the back decay of the muonic-molecular complex [{ddfi)dee]. In this process, 
they gain appreciable kinetic energy. 

Ortho-para effects in /iCF in sohd, liquid, and cold-gas D2 targets with different con- 
centrations of ortho-D2 molecules have been experimentally studied pTlfTS] . The measured 
ortho-para effects in solid deuterium targets (stronger resonant ddfi formation in para- 
rich D2) are opposite to those predicted by theory [15], which requires further experimental 
and theoretical investigations. On the other hand, theory and experiment give the same 
sign of the ortho-para effect in the case of cold D2 gases. 
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Figure 6: Time spectrum of neutrons from d-d fusion in ortho-D2 and para-D2 gas at 
T = 40 K and density = 0.2. 



The calculated fusion-neutron time spectra for ortho-D2 and para-D2 gas at T = 40 K 
and ip = 0.2 are shown in Fig. O The two-Maxwell distribution [3j of the initial dfi 
energy has been assumed. After short time, the number of neutrons in ortho-D2 case 
is significantly greater than in the para-D2 target. This difference disappears at much 
larger times. Such a bump in the ortho-rich D2 target at similar conditions has been 
experimentally observed \XS\- Thus, the theory agrees with the experimental findings, at 
least qualitatively. 

In Figs. [714T0] are presented our theoretical results for the D2 perfect-gas target at 300 K. 
The effective resonant formation rates for F = 1/2 and 3/2 are plotted in Fig. [3 These 
functions are more flat than those in Fig. [2] since the thermal distribution of D2 energies 
is much broader at 300 K and more initial rotational states are significantly populated at 
this temperature. Figure [H] shows the Monte Carlo time evolution of the population of the 
diJ,{F = 3/2) atoms at 300 K and = 0.5. The initial statistical value of the upper-spin 
population is equal to 2/3. After about 200 ns, the thermal equilibrium between the lower 
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Figure 7: The rate of resonant ddix formation in d^ scattering from 300-K D2 gas versus 
dn energy. The sohd hues correspond to the upper spin state F = 3/2 of dfi, the dashed 
hues describe the F = 1/2 state. 

and upper spin states is reached. The theoretical downwards spin-flip rate has been scaled 
by the constant factor of 0.6, determined by comparison with the experimental results [20] . 
The thermalization of d/i atoms at such a density is very fast — it takes only a few ns. 
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Figure 8: Time evolution of the spin state F = 3/2 fusion in D2 gas (T = 300 K, density 
= 0.5). 
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Therefore, the initial dfi-eneigy distribution does not affect the fusion-product spectra at 
significantly larger times. 

The calculated neutron time spectrum for this D2 target is plotted in Fig. [9l for very 
large time interval. Apart from the prompt peak, the spectrum can be well described by 
the steady-state kinetics. The Monte Carlo spectra of fusion neutrons and electrons from 
muon decays are shown for short times in Fig. [101 The electron spectrum is flat since the 
muon decay does not depend on the d^ atom energy. 

In conclusion, the theoretical results of Monte Carlo simulations of the /iCF cycle pro- 
cesses in gaseous D2 targets, corresponding to the JINR target conditions, have been ob- 
tained. The theoretical resonant and nonresonant ddfi formation rates and the differential 
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Figure 9: Time spectrum of neutrons from d-d fusion in D2 gas at T = 300 K and at density 
(/) = 0.5. 
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Figure 10: Time spectra of fusion neutrons and of electrons from muon decay in D2 gas at 
T = 300 K and at density = 0.5. 

cross sections for dfi scattering from D2 molecules established the input for our calculations. 
In particular, the thermalization of dn atoms for various targets and the role of nonresonant 
ddfi formation in short-time neutron spectra have been investigated. The nonresonant for- 
mation can be directly observed at very short times in low-temperature targets at densities 
<f) < 0.05. Also, the d/x-thermalization time is important at relatively low target densities. 
The calculated neutron time spectra for ortho-D2 and para-D2 at 40 K and (p = 0.2 dis- 
play the enhanced neutron emission in the ortho-rich targets at intermediate times, which 
agrees qualitatively with the recent experimental data. The simulation for the high density 
(0 = 0.5) D2 target at 300 K shows that the neutron time spectrum can be well described 
by the steady-state kinetics, apart from very short times. 
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